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1 Abstract 

A theoretical approach with a microscopic model is proposed for the observed "high 
temperature superconductivity" in the Iron-based compounds. The above scheme takes 
into account two important aspects viz. (i) superconducting transition close to magnetic 
ordering and (ii) the layered structure. From the calculation of the superconducting 
transition temperature, it is shown that in the Fe-based superconductors the magnetic 
mechanism for superconductivity, operating through the effective attractive Coulomb in- 
teraction within the framework of the Fermi Liquid theory, is highly plausible. 
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2 Introduction 



There has been a lot of excitement over the recent discovery of superconductivity in a 
new class of systems viz. the doped phase of layered rare earth-transition metal based 
compounds in the temperature regime of 25-50K. They have a lot of similarities with the 
high temperature superconductors belonging to the Cuprates family in certain important 
experimental features. These include (i) layered structure with specific roles of various 
layers and (ii) the proximity to a magnetically ordered state. In this article we will only 
focus on the systems which contain Fe as the transition metal component. 

These Fe-based superconductors can be divided into 3 types depending on their chem- 
ical composition. They are :- (1) RO\- x F x FeAs where R represents a rare earth element 
like La, Sm etc.; (2) X\- x Y x Fe2As2 where X stands for any bivalent metal like Ca, Ba, 
Sr etc. and Y represents a monovalent metal like Na, K etc. and (3) XFe\~ x Y x AsF 
where X implies a bivalent metal and Y is a transition metal different from Fe. The 
experiments have shown that in all these systems the Fe-As layer is responsible for both 
electrical transport as well as magnetism and that the magnetism is of itinerant nature 
[1,2]. In the parent or undoped phases, i.e. the phases corresponding to x — 0, all the 
above systems are insulator and exhibit long range magnetic ordering of spin density wave 
type. For the systems belonging to the subclass (1) above , the rare earth-oxygen layer 
acts as a charge reservoir and provides electrons to the Fe-As layer to act as carriers when 
O is partially replaced by F. This turns the system conducting and also suppresses the 
long range magnetic ordering. The system behaves as a normal paramagnetic metal with 
enhanced spin correlations and becomes superconducting at lower temperatures. The su- 
perconductors belonging to subclass (2) are hole doped and those belonging to subclass 
(3) can be either electron or hole doped. The common feature exhibited by all these 3 
subclasses of systems is however emergence of superconductivity close to magnetic order- 
ing under doping. 

It is worthwhile to point out that the behaviour of the normal conducting phases of 
these Fe-based systems is very close to that of a usual metal, unlike those of the several 
members of the Cuprate family. Moreover, the complete phase diagram of these Fe-based 
systems as a function of doping, is much less complex compared to those of the Cuprates. 
Besides, from the phenomenology it seems that the roles of phonons and other electronic 
charge excitations on superconductivity is rather unimportant [3]. 

On the basis of these above observations, we propose a simple microscopic model here 
for the superconducting phase of Fe-based superconductors. In our theoretical model we 
invoke a mechanism for electron pairing close to magnetic ordering, based on the dielec- 
tric function formalism [4]. Moreover, we also incorporate the crucial role played by the 
anisotropic layered structure in our calculations to some degree. We then explore the 
possibility of the effective Coulomb interaction between the electrons (or holes) becoming 
attractive due to the presence of the enhanced spin fluctuations in the system. Finally, 
we present an analytical formalism for calculating the superconducting transition tem- 
perature for our system, when the above mechanism is operative. 
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3 Mathematical Formulation and Calculation of Su- 
perconducting Transition Temperature 

We explore the pairing scheme of Bardeen-Cooper-Schrieffer (BCS) in a Fermi Liquid 
(FL) background. The assumption of the FL description of the normal phase is justified, 
as r s values of the equivalent electron (or hole) gas in Fe — As layers in the doped phases 
of these systems fall in the range of 4 to 5, as determined from the experimental values 
of the carrier concentration and furthermore there is no experimental support for strong 
on-site Coulomb correlation either [5,6]. 

(A) We first propose a very generalized model which can be appropriate for these 
Fe-based superconductors. This is written as :- 

Kgen. = Y e KrnC£a,m C ka,m + Y ^kk',m&k',rrA,m + Y ^k' [&k',nA,n + h.C.] (1) 
m,k k,k',rrt <mra>,k,k' 

where, in the right hand side the first and the second term together represent the standard 
BCS Hamiltonian corresponding to the intra-layer pairing with m being the layer index 
for a particular Fe-As layer. The e(k) dispersion is highly anisotropic with the major 
contribution coming from the intra-layer kinetic energy. The c's are the usual single 
fermion operators and the 6's are the Cooper pair operators. The last term represents the 
hopping of Cooper pairs between two successive Fe-As layers. It may be recalled that fe's 
are related to the c's in the following way, 

K = c_ k _ CT c k(7 (2) 

The basic idea is that the Cooper pairs are formed in each of the above mentioned layers 
through an effective attractive interaction Vkk',m ; however since superconductivity is a 
3-dimensional phenomenon, we had to invoke an inter-layer process coupling the pairs 
from different layers with a parameter A. This inter-layer process may also contain a 
screened Coulomb interaction due to the background dielectric constant (with contribu- 
tions from the layers containing the rare earth atoms), accompanying the pair tunnelling. 
It should be stressed that the parameter A is expected to be much smaller compared to 
V. Furthermore, we assume that the intra-layer attractive pairing interaction is a func- 
tion of the in-plane momentum transfer of the electrons i.e. k|| — kj| only. The pairing 
interaction can be determined from the longitudinal dielectric function of the 2d electron 
(or hole) gas, which can take into account the spin fluctuations. Besides, the out of plane 
components of the electron momenta are taken to be box quantized in the following form 

k - ^ (3) 

^perp. \ U J 

where p is an integer, +ve or -ve including zero and L is the size of the lattice in the z 
direction. 

We could define the order parameters viz. the BCS gap function and the pair tunnelling 
function in the following way 

Ak = -E^kVn<#', m > ( 4 ) 
k',m 
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and 

T k = E AS? < 6£ >m > (5) 

k',m,n 

Carrying out a standard mean field treatment on our Hamiltonian, we generate the 
coupled self-consistent equations for the above order parameters. Our preliminary calcu- 
lations with these bring out interesting results for the Bogoliubov quasi-particle excitation 
spectra and the superconducting gap equation. The quasi-particle spectra determined for 
our anisotropic layered system, exhibits significant and non-trivial departure from the 
conventional one known for the isotropic 3d system [7]. In particular, the expression for 
the excitation energy gap turns out to be slightly different from the BCS gap parameter. 
Detailed calculations are in progress and the results will be reported later. 



(B) We now work with a simpler version of the earlier model i.e. the generalized model 
given by equation (1). Here we consider only the first two terms from equation (1) and 
remove the explicit layer indices. Therefore our simplified model is similar to the BCS 
model for an isotropic system; however the pairing interaction V is once again assumed to 
be a function only of the in-plane momentum transfer and is explicitly calculated from the 
response functions corresponding to an interacting 2d electron (or hole) gas. In particular, 
this dielectric function includes exchange-correlation corrections beyond Random Phase 
Approximation (RPA) and can be related to the spin susceptibility [4,8]. The repulsive 
bare Coulomb interaction Vb(q) is also taken corresponding to the 2d electron gas and is 
given by 

^o(q) = ^ (6) 

where, e# is the background dielectric constant of the system, which has the main contri- 
butions from the layers containing the rare earth ions, besides Oxygen and Flourine ions 
for some of the Fe-based superconductors. 

Let us focus on the relevant set of equations now. From the expressions for the static 
and wave-vector dependent charge and spin responses for 2d electron (hole) gas, we arrive 
at the following important equations [4,8] 

e-\q) = 1 - V (q)n(q); n'^q) = n^(q) + V (q)[l - G s (q)\ (?) 

where e(q), ir(q) and 7r (g) are the longitudinal dielectric function, the full polarizability 
function and the RPA polarizability function respectively. The function G s (q) is the spin 
symmetric local field correction. The expression for 7r (g) is given by [9] 

MQ) = ^ptk'l + 2sgn(u)9(u 2 - \){v 2 - 1)*] (8) 

where |g'| = v = k F is the Fermi wavevector and N(0) is the single fermionic 

(electron or hole) density of states for one kind of spin at the Fermi surface. Similarly, 
we have the following expression for the magnetic spin susceptibility 
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1 - V (q)G a (q)n (q) 
where G a (q) is the spin anti-symmetric local field correction. 



From the above set of equations (7) - (9) we arrive at the following relation connecting 
the electrical and the magnetic response [4,8] 



The above equation (10) brings out an extremely important quantum effect. It shows 
that for the q values where the magnetic spin susceptibility is highly enhanced over the 
magnitude expected from the RPA or Pauli response so that the first quantity (in the 
bracket) on the right hand side becomes negative, then the dielectric function corre- 
sponding to those q modes can itself become negative. It can be argued that the function 
7r _1 (g) is always positive because of the stability criterion [4]. This has the immediate 
consequence for the effective static Coulomb interaction given by V e ff(q) = e~ 1 (q)Vo(q) 
viz. the effective interaction becoming attractive corresponding to the q modes where a 
paramagnetic layer exhibits enhanced magnetic spin susceptibility. This further leads to a 
scenario where the Cooper pairing can take place through this electronic mechanism itself. 

Let us now relook at the experimental situation for the Fe-based superconductors. 
The parent systems order magnetically with a long range SDW with q of the order of kp- 
In the doped phase which is paramagnetic and metallic, the magnetic spin susceptibility 
is highly enhanced in the above q regime, particularly when the doping level is not too 
high. The results from neutron scattering experiments and ac susceptibility measurements 
confirm this [10]. Therefore, our proposed electronic mechanism for superconductivity in 
the Fe-based systems is indeed very very realistic. 

We would now like to explore the feasibility of such a mechanism to yield the super- 
conducting transition temperature (T c ) in the range as observed for these superconductors. 

In the weak coupling BCS theory for superconductivity in the s-wave channel, the 
attractive coupling constant is obtained by averaging the static g-dependent attractive 
interaction over the Fermi surface [4,11]. The equation for T c is related to this coupling 
constant and the cutoff energy (temperature) in the BCS square well model for the ap- 
propriate pairing interaction. In our case, this interaction is of the attractive Coulombic 
type and thus the cutoff scale is the Fermi energy (temperature). 

The attractive coupling constant (A e ) in our case is given by 



= \Ax-\q) ~ V (q)(G s (q) - G a (q))]7T~ 1 (q) 



(10) 




(11) 



provided the right hand side of the above equation is negative; kF\\ is the Fermi wavevec- 
tor corresponding to the effective 2d electron gas, which is related to the magnitude of 
the corresponding r s , referred to earlier. 



The equation for T c under the above mentioned mechanism, then becomes 

T c = l.UT F exp{^-) (12) 
|A e | 

where Tp is the Fermi temperature of the system. 

Assuming the value of E F to be between 0.6 ev to 1.0 ev for these Fe-based systems in 
the metallic phases, as determined from the photoemission experiments [6], we get from 
the above equation (12) the theoretical range of the magnitude of the attractive Coulomb 
coupling constant lying between —0.1 to —0.2 corresponding to the regime of observed T c 
in the broad spectrum of 25K to 50i^. Thus it is quite possible to attain the transition 
temperatures in this class of superconductors by invoking the magnetic mechanism in the 
weak coupling regime. 

We now examine more carefully the realisation of this strength of the attractive cou- 
pling from the microscopic many-body parameters used in our analysis earlier. Detailed 
calculation involving equations (6)- (11) and making use of the results from the earlier 
works [4,8,9], it can be shown that an average magnetic spin susceptibility enhancement 
over the g-space of the order of 10 or even slightly less, can easily enable us to reach 
the above magnitude of the attractive Coulomb coupling constant. In this study, the 
typical bare repulsive Coulomb coupling constant (obtained at the RPA level, neglecting 
the local field corrections) was taken around +0.3 and the g-space averaged ratio of the 
very important and crucial correlation parameter Vo(q)[G s (q) — G a (q)} and the RPA spin 
susceptibility, was kept within the allowed range for the metallic density limit. This the- 
oretically obtained magnitude for the averaged spin susceptibility enhancement is very 
plausible, considering the fact that the strongest enhancement in these experimental sys- 
tems takes place for the modes with q around k F , as the parent systems exhibit the long 
range SDW ordering. 



4 Conclusions 

In this article, we have tried to present a simple but realistic microscopic model and 
approach for investigating the superconductivity in the new class of Fe-based super- 
conductors. We have proposed a magnetic mechanism for this layered system and have 
considered the conventional s-wave pairing. An approach very similar to the present one 
was applied by us successsfully to a ferromagnetic superconductor I9C07 earlier [4]. The 
scheme presented here, can easily be extended for exploring the anisotropic pairing chan- 
nels within the BCS or Eliashberg scheme as well. 

We have made use of the many-body parameters extracted from the electron gas cal- 
culations. More accurate quantitative prediction can only be made by taking into account 
the detailed band structure effects. 

Many of the member systems in this family exhibit a structural transition in the 
vicinity of the long range magnetic ordering in the undoped parental phase [1-3]. This 
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gives rise to the possibility of the availability of soft phonons even in the doped metallic 
phases. These phonons may also play a secondary role in mediating the pairing interac- 
tion, besides the major contribution coming from the magnetic mechanism described here. 

Our theoretical approach can be strengthened by more detailed experiments. In par- 
ticular, a study of the frequency integrated spectrum of the dynamical structure factor 
S(q, u) in the doped phases, extracted from the inelastic neutron scattering experiments 
can throw lot of valuable information regarding the g-space averaged spin susceptibility 
enhancement. 
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